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Using a highly efficient Monte Carlo algorithm, we are able 
to study the growth of coverage in a random sequential ad- 
sorption (RSA) of self-avoiding walk (SAW) chains for up to 
~ 10 12 time steps on a square lattice. For the first time, the 
true jamming coverage (8j) is found to decay with the chain 
length (N) with a power-law 8j oc AT -0 ' 1 . The growth of the 
coverage to its jamming limit can be described by a power-law, 
8(t) ~ 8 j — c/t v with an effective exponent y which depends 
on the chain length, i.e., y ~ 0.50 for N = 4 to y ~ 0.07 for 
N = 30 with y — > in the asymptotic limit N —* oo. 
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Studying the kinetics of random sequential adsorption 
(RSA) has attracted a considerable interest in recent 
years [flj— pTsfj because of its enormous applications 1 19 24|] 
in the adsorption processes involving a variety of species 
from a point-like particle to a protein-like complex struc- 
ture in physical, chemical, and biological systems. Some 
of the examples include binding of ligands on polymer 
chains, coating, designing composites, chemisorption, ph- 
ysisorption, and reaction of molecular species including 
globular protein on surfaces and interfaces, etc. These 
adsorption processes may be divided into two categories: 
(1) annealed adsorption where the species are mobile (a 
thermal equilibration for the interacting adsorbants) be- 
fore they settle onto the surface - a cooperative sequential 
adsorption. (2) Quenched adsorption where the adsorp- 
tion occurs without subsequent diffusion or desorption. 
We consider the latter category known as random se- 
quential adsorption. 

The problem of RSA in one dimension p9|j25| -p7| is 
well understood with exact results for some adsorption 
processes. Understanding the growth of coverage in two 
dimensions with the RSA lacks rigorous results by ana- 
lytical methods due to their intractabilities especially for 
objects with polydisperse shapes. Therefore, computer 
simulations remain one of the primary tools to investi- 
gate these problems. Numerical results || and theoreti- 
cal analyses for the deposition of disks on continuum 
suggest that the coverage follows the Feder's law 0] at 
large time, 



0(t) « Oj 



c 



(1) 



where, 9(t) is the coverage at time t, 9j is the jamming 
coverage, d is the dimensionality of the host space, and 



c is a constant. 

The temporal dependence, Eq. ([!]), seems to describe 
the deposition of hyperspheres. For squares H, rectan- 
gles M, and ellipses Q, on the other hand, the above 
law is still valid but d is replaced by d/, the number of de- 
grees of freedom of the corresponding objects. For disks 
with polydisperse sizes the computer simulation studies 
|l3| seem consistent with the theoretical prediction jl2| , 
9(t) » 9 j — ct~ v , where y = l/(d+ 1). As a natural ex- 
tension, one would like to consider objects not only with 
polydisperse sizes but also with polydisperse shapes, a 
more complex problem in the studies of RSA. Therefore, 
we perform a large-scale Monte Carlo simulation to study 
the adsorption of polymer chains which have polydisper- 
sity in both shapes and sizes. A model chain is one of the 
most simple ramified objects with its well-known shape 
and size distributions in a variety of systems p8pfj| ] . The 
RSA of polymer chains is nevertheless relevant in applica- 
tions such as coating and paint. In contrast to theoretical 
predictions for the RSA of polydisperse objects of regular 
shapes, we find that an effective exponent y depends on 
the chain length. Furthermore, we are able to reach the 
true jamming coverage (9j) with an efficient algorithm 
which enables us to predict a power-law dependence on 
the chain length (N), 9j oc JV -0,1 . 

We consider an L x L square lattice with a periodic 
boundary condition. A polymer chain is modeled by a 
self-avoiding random walk (SAW) which is generated on 
the trail of a non-reversal random walk (NRRW) with 
self- avoiding constraints. The chains are dropped onto 
the square lattice, one at a time sequentially. If a chain 
overlaps with previously deposited chains, the attempt is 
rejected. Once the chain is deposited on the lattice, it 
sticks on the surface permanently. The deposition rate 
or equivalently the time scale is fixed with the follow- 
ing algorithm. In unit time At = 1, L attempts are 
made to generate and deposit SAW chains each of length 
N starting from a randomly selected site. As soon as 
a walk overlaps with the previously deposited site (from 
other chains or itself), the walk is abandoned, and a new 
attempt is made to deposit a chain starting from a ran- 
domly chosen new site. 

Due to the large number of conformational states of 
polymer chains, the deposition becomes very slow in the 
late stage. An event-driven method |H]|j(]|3lJ is, there- 
fore, used to speed up the simulation. This is accom- 
plished by identifying the early part of growing chains, 
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i.e., the partial chains and then classifying them peri- 
odically as "available" or "not available" for the deposi- 
tion in an iterative fashion. The subsequent depositions 
are then made starting from the available partial chains 
choosing at random. The list of the partial chains is lim- 
ited by the available computer memory and governs the 
speed of the program. It appears that there is an optimal 
list size for a given value of N. 

Now, in order to maintain the same dynamics, each 
deposition attempt does not simply advance the time by 
St = L~ 2 as before, but by an amount jSt, where j is a 
random integer with an exponential distribution, p(l — 
pY -1 . p is the ratio of the number of potentially available 
partial chains in the list to the total number of chains of 
partial length n, L 2 Z™ rrw , where Z™ rrw is the number of 
NRRW chains of length n. The random integer can be 
generated by 



J 



1 



ln(l-p) 



(2) 



where £ is a uniformly distributed random number be- 
tween and 1. 

The simulation is performed on IBM SP2 and fast 
workstation clusters. The total amount of CPU time 
for the computation is about six months equivalent of a 
single DEC AlphaStation 250/266. Variation of the cov- 
erage 9 with time t is presented in Fig. [l]for various chain 
lengths. We immediately notice that a rapid increase in 
the coverage in the short time regime is followed by a very 
slow growth in the long time. We can divide the char- 
acteristic behavior in three time regimes: (1) short time 
regime, (2) intermediate regime, and (3) very late-stage 
regime. 

The short and very long time regimes (1) and (3) are 
easily understood. The coverage at very short time is pro- 
portional to t, since the depositions of nearly all SAWs 
are accepted. In very late stage [ |32| , the dynamics is con- 
trolled by filling independently the last pores. The holes 
disappear according to e~*/ T , where 1/r is the probabil- 
ity that a given hole is being filled (in unit time). For the 
deposition of chains, the slowest mode is filling the void 
in exactly two ways (out of Z^f rw possibilities). Thus 



*ynrrw 
N 



2-3 



N-2 



N > 1. 



(3) 



The exponential decays and the prediction for r agree 
very well with simulation data, see Fig. |^. 

In the intermediate time regime (2), 1 < t <C t, a 
power-law dependence is observed even on discrete lattice 



i.e., 



9{t) 



ct~ 



(4) 



with an effective exponent y. Since the jamming limit 6j 
is not known accurately for long chains, it is instructive 
to consider the derivative of 8 with respect to t, 



(W _L_ 

dt K ti+i/" 



(5) 



This quantity is plotted in Fig. ^| for various chain 
lengths. We see a crossover from a power-law variation 
of the rate of coverage in the intermediate time regime 
to an exponential decay in the long time especially for 
short chains. Note that this intermediate regime expands 
on increasing the chain length. We observe a remark- 
able power-law behavior over 12 decades in time with 
the largest chain length N = 30 for which the true jam- 
ming limit is reached in our simulations. Least-square 
fits in the power-law regime give a reliable estimate of 
the effective exponent y, particularly for N < 30. Note 
that the small value of y ~ 0.07 for N = 30 is stable 
over a relatively large time scale within the statistical 
fluctuations (see the inset in Fig. ||, and Fig. ||). A loga- 
rithmic fit 8(t) » 6 j — c/log(t) appears less satisfactory. 
We see that the magnitude of the exponent y depends 
systematically on the chain length, y ~ 0.50 for N = 4 
to y ~ 0.07 for N — 30. To our knowledge, none of the 
previous studies has shown a size-dependent exponent y. 
A crude extrapolation (see Fig. |^) leads to a small value 
of y for larger chain length with a possibility for y — > 
as N — ► oo within the statistical error. 

Let us recall the well-known Swendsen's argument 0] 
for Feder's law [Eq. (|j])] for the adsorption of disks on 
continuum. At late stage, the pore vanishes according 
to e~ kt , where k oc l d , I is the linear size of the pore, 
and d is the spatial dimension. Assuming an uniform 
linear size I, the probability distribution for the rate k 
is p(k)dk oc 1 • dl oc k x l d ~ x dk. Total contribution to the 
approach to jamming is then 



0{t) oc / k^- 1 e- kt dko^r 1/d . 



(6) 



Viot and Tarjus j7j generalized the above result to 
anisotropic objects and concluded that the approach 
to jamming is t~ 1 / d f for monodisperse objects, where 
df is the degrees of freedom of the objects. For the 
objects with regular shape but polydisperse sizes p2| , 
9;-J(t)ar 1 /! J /+ 1 ), 

Polymer chains have much more internal degrees of 
freedom. For sufficiently long chains, the lattice structure 
becomes less important. The chains can be specified by 
the orientation of their segments on a coarse-grain level. 
The number of degrees of freedom is proportional to the 
number of such segments. Thus, we may identify df oc N 
and consequently y oc 1/N. In fact the data for y can be 
roughly characterized by y ss 2/N. 

Alternatively, one may identify the degree of freedom 
as the number of SAW configurations Zf^ w of length N, 



i.e., df oc Z 



ATT- 



-1~N 
Z eff 



where 7 is a critical expo- 



nent and z e ff < z — 1, where z is the coordination num- 
ber @. Then, the exponent y « 1/{Z S ^ W + 1) becomes 
much smaller than our estimates. Thus, the application 
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of Tarjus and Talbot result JT|] is not valid if d f = Zff w 
is assumed. However, in asymptotic limit for large N, 
y — > as N — > oo is consistent with the prediction of 
Tarjus and Talbot. 

With the event-driven method we are able to reach 
jamming much faster which gives us accurate estimates 
for the jamming coverage Oj, see Table | The jamming 
coverage decreases with the chain length. It is not clear 
a priori that the jamming coverage 6j(N) goes to zero as 
N approaches infinity. However, assuming. 

6j(N) cx N~ x , (7) 

a Oj(N) versus N plot on a log-log scale (see Fig. ||) 
leads to a very small exponent x w 0.11. Note that the 
magnitude of the exponent x is much smaller than that of 
ref. |[6| (x « 1/3) which deals with the percolation f34f| 
of chains. The jamming coverage in the percolation, i.e., 
the maximum value of the volume fraction is the "first 
jamming coverage" which is far from the asymptotic limit 
of our "true jamming coverage" here. 

We have also analyzed the conformation of the chains 
during the deposition process. For SAW, it is well-known 
p8| , p9| that the radius of gyration or end-to-end distance 
varies with chain length as 

R g cx R e cx N 3/4 . (8) 

The question now is whether this law is modified for the 
chains deposited on the lattice. Our numerical data in- 
dicate that the average radius of gyration is less sensitive 
to packing. One of the possible reasons for not observing 
much change in the power-law dependence from that of a 
free SAW, is the fact that the conformations of the chains 
in their early stage of deposition dominate the average ra- 
dius of gyration. However, due to the irreversible nature 
of the problem, the chains have memory of their arrival 
time in the sense that their conformation and correlation 
depend on t. In fact, we do observe that the radius of 
gyration decreases systematically with time especially in 
the long time regime. 
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FIG. 1. Coverage as a function of time, for chain length TV 
indicated by the number. See table Q for the statistics on the 
sample size and the number of independent runs. 

FIG. 2. Semi-logarithmic plot for the difference between 
jamming coverage and coverage at time t, for chain lengths 
TV = 2, 5, 10, and 15. For a better view of all the curves 
on the same plot, we plotted again the normalized time t/r. 
The straight line is the pure exponential decay e _t for TV = 1. 
The insert shows r as a function of TV, the straight line is the 
theoretical prediction, Eq. (||). 

FIG. 3. Derivative of coverage vs. time on a log-log scale. 
The number indicates chain length TV. Insert shows the vari- 
ation of exponent y with time evaluated at a regular interval 
for TV = 30. The statistics is presented in table ffl. 

FIG. 4. Exponent y obtained by least-squares fits of Fig. ^| 
plotted against inverse chain length. 

FIG. 5. Jamming coverage versus chain length on a loga- 
rithmic scale. 



TABLE I. Jamming coverage 6j = 9(t — > oo) for various 
chain lengths. Statistical errors on the last digits are indicated 
by the numbers in the parentheses. 



TV 


$ j (error) 


L 


number of runs 


1 
2 


1 

0.906820(2) 


1000 


26200 


3 


0.858296(4) 


500 


34000 


4 


0.837055(13) 


500 


14000 


5 


0.81235(1) 


500 


7501 


7 


0.78558(2) 


500 


2320 


10 


0.75895(1) 


500 


4701 


15 


0.72473(8) 


500 


700 


20 


0.70178(10) 


500 


33 


30 


0.6683(4) 


500 


28 


40 


0.654(9) 


200 


7 
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Fig. 2 




